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Abstract. We propose an experiment to measure the electric dipole moment of the 
electron using ultracold YbF molecules. The molecules are produced as a thermal 
beam by a cryogenic buffer gas source, and brought to rest in an optical molasses that 
cools them to the Doppler limit or below. The molecular cloud is then thrown upward 
to form a fountain in which the electric dipole moment of the electron is measured. 
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1. Introduction 

The standard model of elementary particle physics predicts that the electric dipole 
moment (EDM) of the electron is exceedingly small fT], c::^ 10~^^e.cm, as illustrated 
in figure [T] In all local, Lorentz-invariant field theories a permanent EDM requires CP 
symmetry to be violated. This happens in the standard model through the Yukawa 
couplings in the quark sector [2j, with the first non- vanishing result appearing at the 
four-loop level [3j. However, the standard model is thought to be incomplete, one 
reason being that it does not have enough CP violation to explain the excess of matter 
over antimatter in our universe, as first discussed by Sakharov [4]. It would seem that 
new forces are waiting to be discovered, involving new particles with masses probably 
near the electroweak scale of 0.2 TeV up to a few TeV. This motivates many current 
experiments in particle physics. 

The electron is sensitive to these new CP- violating forces because they are expected 
to induce an EDM, even at the level of a simple one-loop radiative correction, such as 
the one shown in the upper right of figure [!} Here the natural size of the EDM is [5j 
sin(0(7p), where a is the fine structure constant, mg is the mass of the electron, 
Ms is the mass of the new particle, here the selectron, and 9c p is the CP violating 
phase of the coupling. Since there is no reason to assume CP conservation, let us take 
Ocp , then M, = ITeV/c^ gives ^ 2 x IQ-^^e.cm. Similarly, = lOOTeV/c^ 
gives de — 2 X 10~^^e.cm. In this example the virtual selectron and gaugino are 
supersymmetric particles, but most other models similarly predict EDM values ^ in 
the range 10~^^ — 10~^^e.cm, as shown in figure in If there is new particle physics in the 
energy range 200 GeV-100 TeV, an electron EDM sensitivity of 10~^°e.cm is very likely 
to detect it. 

Experiments going back to the 1960s [8] have sought to detect the electron EDM 
using beams of heavy atoms. These have the virtue that the interaction energy rjdeE 
between the electron EDM d^ and the applied electric field E is enhanced by the 
relativistic factor ry, which can be large, as discovered by SchiflF [9] and Sandars [lOj. 
The atomic beam experiments ended in 2002 with a measurement on Tl, a complex 
tour de force that gave the limit [Hj \de\ < 1.6 x 10~^''e.cm. In order to improve on 
this, a new technique has been developed using polar molecules [12], which oflFer two 
great advantages, (i) Molecules are more sensitive to de than atoms. For example, 
in a large laboratory field YbF molecules are 200 times more sensitive than Tl atoms 
[13] . (ii) Molecules can be immune to the particular systematic error (generated by the 
motional magnetic field E x v/c^ ) that limited the Tl experiment [H]. The recent 
electron EDM measurement using this new method [15] has given a slightly improved 
upper limit of \de\ < 1.05 x 10~^''e.cm, shown by the dashed line in figure [!} Together, 
these measurements indicate that the simple supersymmetric model illustrated in figure 
[T] cannot be right: either the CP violation is suppressed for some reason or the mass of 
the new particle is considerably greater than 1 TeV. While it is possible to construct 
supersymmetric theories where some particular EDM is accidentally small, it does not 
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Figure 1. Coloured boxes: range of predicted values for the EDM. Dashed line: 
current upper limit on electron EDM [15]. Solid line: sensitivity of proposed new 
method. Top right diagram: example of a simple supersymmetry 1-loop radiative 
correction, generating an electron EDM. A cross marks the CP- violating vertex. 



seem possible in a minimal supersymmetric model to keep the electron and neutron 
EDM both small enough [16] . One is led to the conclusion that either the superpartner 
masses are rather heavy or there is some more complex physics on the electroweak scale 
that breaks supersymmetry but manages to suppress CP- violation. This is an important 
open question at a forefront of fundamental physics that will be elucidated by a more 
sensitive search for the electron EDM. 

To go significantly beyond the current level of EDM sensitivity, one needs a 
fundamental improvement in the method, which is what we propose here. By laser- 
cooling the molecules, it will be possible to replace the molecular beam by a molecular 
fountain, where each molecule can precess coherently for almost one second, rather 
than the current one millisecond. This will allow detection of an EDM as small as 
1 X 10~^^e.cm, corresponding to CP-violating physics at energies up to 100 TeV. With 
this sensitivity one can surely hope to uncover some of the new physics underlying the 
important issue of CP-violation in the early universe. 

Quite apart from this application in elementary particle physics, there is strong 
interest in cooling molecules to microKelvin temperatures for possible use in quantum 
information processing [TTl [IE], simulation of quantum many-body systems [191 ED], 
ultracold chemistry [211 ES] , nuclear physics ^23j , and cosmology [24j . 

One approach is to photo-associate or magneto-associate atoms that are first 
separately cooled to ultra-low temperatures. This is proving to be a powerful tool 
for studying the physics of ultracold molecular gases [25l [26]. However, the method 
is limited to a few alkali/alkali and alkali/alkaline-earth diatomics and is therefore not 
suitable for many of the applications listed above that require other types of molecule. 
By contrast, a huge range of molecular species can be produced by supersonic expansion 
or by cryogenic buffer gas cooling [i27j . These methods produce molecular beams whose 
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Figure 2. Scheme for laser cooling YbF. The cooling is due to multiple scattering of 
552 nm light on the - X^I]+(v = — 0) transition. Molecules that decay iov = l 

and V = 2 are quickly repumped to v = using lasers at 568 nm and 584 nm. Molecules 
that fall into the v = 3 state are cold and require only weak repumping. The branching 
ratios to the various vibrational levels, and the lifetimes for spontaneous vibrational 
decay in the X state, are also shown. 

temperatures are typically a few Kelvin. This is not cold enough for the proposed 
molecular fountain which requires a temperature in the microKelvin range. For this 
reason, it is important to find a cooling method that bridges the gap between a few K 
and 100 /xK. 

One possibility for bridging this gap is to trap the molecules together with ultracold 
atoms and allow them to cool sympathetically as a result of their collisions with the 
atoms pHl EHl [30l EU |32j. It is possible that this method will yield very cold, dense 
molecular samples if we can find cases where the elastic collision cross section is large 
enough and the losses are small enough. However, the number of molecules may well 
be inadequate for EDM measurements, which require low statistical noise and must 
therefore detect many molecules. 

Here, we propose the promising alternative of using laser cooling, a method already 
well established for atoms. Most molecules cannot be laser-cooled because they decay 
after a few spontaneous emissions into a higher vibrational state that does not interact 
with the cooling light. However, there is a family of molecules where this difficulty is 
not too severe ^33j. De Mille's group at Yale has recently demonstrated laser cooling of 
SrF molecules [SU [35]. In the case of YbF, our group has measured the Franck-Condon 
factors [36j, the ratios that characterise the change of vibrational state, and find them 
to be suitable for laser cooling. 

Figure [2] shows the relevant energy levels of YbF. The lowest odd-parity state 
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Figure 3. The main components of the proposed experiment. 



X^i;+(^ = 0,A^ = 1) is excited on the P(l) line to the even parity k^Ili/2{v = 0) 
state by near-resonant hght of wavelength Aqo = 552 nm. This can only decay to X^S+ 
states having = 1 because the next odd-parity ground state, = 3, is forbidden by 
the angular momentum selection rules for electric dipole radiation. Each scattered 
photon imparts a momentum /^/Aqo to the molecule {h being Planck's constant), 
corresponding to a velocity change of 3.7 mm/s. By comparison, YbF molecules in 
thermal equilibrium at a temperature of T = 3.6K have an rms velocity along a given 
direction of ^Jk^TJM — 12m/s, where fc^ is the Boltzman constant and M is the mass 
of one molecule. Thus, molecules pre-cooled to this temperature by a helium buffer gas 
can be brought to rest by scattering a few thousand photons. Six beams propagating 
along the Cartesian directions ±x, ±z are detuned by approximately 5 MHz to the 
red of resonance in order to form an optical molasses [37]. This cools the molecules to 
a temperature of order hT/{2kB) = 140/iK, where F is the spontaneous decay rate of 
the A state. There is a 6.6%(0.3%) probability that the electronic decay will excite the 
V = l{v = 2) vibrational state [36j. These molecules are returned to the cooling cycle by 
repump lasers at 568 nm (584 nm). We show below that a 3.6 K molecular beam can be 
cooled to form a molecular fountain and that this will improve the experimental EDM 
sensitivity by a very large factor, as indicated by the solid arrow in figure [T] 

Figure |3] illustrates the main components of the proposed experiment. The YbF 
molecules are first produced as a pulsed 3.6 K beam through thermalisation with a 
cryogenic helium buflFer gas. They are then steered away from the helium beam and 
into a region of lower background gas pressure by a magnetic guide. On exiting from 
the guide each cloud of molecules encounters a region of optical molasses, where it is 
stopped and cooled to ultralow temperature. The cold molecules are then thrown up 
as a fountain into the space between two electric field plates. There the electron spin 
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Figure 4. Cryogenic, thermal molecular beam source design, taken from [38] 



is first polarised, then allowed to precess in the electric field. Finally, the spin direction 
is read out and the electron EDM is deduced from the amount of precession. In the 
following sections we describe each of these steps in more detail. 

2. Thermal beam 

The YbF molecules will be prepared as a thermal beam at a temperature of 3.6 K, using 
the double cryogenic cell design developed by Doyle's group [38] and illustrated in figure 
[4} Helium gas enters a first cryogenic cell, where a target containing Yb metal and AIF3 
is illuminated by pulses of infrared light from a YAG laser. The materials ablated from 
the target react to produce a cloud of YbF molecules. These thermalise with the cold 
helium gas but emerge through the first aperture with a boosted centre-of-mass velocity 
due to the fiow of the helium. The resulting velocity distribution is well described by 
f{u) oc 2i^exp(— M(2i — UqY /{ksT))^ where Uq is typically 100 — 200 m/s. The second 
cell, where the pressure and helium fiow velocity are much lower, reduces Uq close to 
zero. The collisions in this cell reduce the total number of molecules in the beam, but 
greatly increase the number of slow molecules as shown experimentally in |38j. This 
is important because the laser cooling cannot bring faster molecules to rest. In our 
own group, we have studied the production of YbF in a single cryogenic buffer gas cell 
|39lll0]. We know from these experiments that we can produce 1 x 10^^ X-state molecules 
per pulse in the cell [40] at the He density that optimises thermal beam production. We 
have extracted a YbF beam from such a cell [H] , but have not tried making a thermal 
beam using a two-cell source of this type. Reference [38] reports an extraction efficiency 
of 0.5% for CaH, and we see no reason why the same should not apply to YbF. With 
this assumption, we expect a thermal beam of 5 x 10^° X- state molecules per pulse at 
a temperature of 3.6 K, of which 24% (according to the Boltzmann distribution) will be 
in the state X'^Ti^{v = 0, = 1) that is to be laser cooled. These numbers are recorded 
for later reference in Tabled! 
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A 20 cm long magnetic guide collects the molecules that emerge from this source and 
bends them out of the helium beam. The bend radius is large enough that only a 
small fraction of the very fastest molecules are lost, and these cannot be captured by 
the molasses anyway. Magnetic guiding from a cryogenic effusive beam has previously 
been demonstrated using O2 molecules [42j. We propose the same type of guide: a 
magnetic octupole built from Nd:Fe:B magnets. The eight poles are arranged with their 
faces touching a circle of radius 5 mm, and with gaps of 1mm between their corners. 
According to a numerical model of the magnetic field, this trap has a depth of 0.6 T 
and therefore confines weak-field-seeking molecules whose transverse velocity is below 
5.9 m/s, amounting to 13% of the Boltzmann distribution. Since half the molecules are 
in weak-field-seeking states, the transmission of the guide is expected to be 6.5%. There 
should be no significant loss of molecules from collisions with the helium because the 
cloud coming out of the cell is thermal, indicating that the molecules are decoupled from 
the much faster helium beam. Spin fiip losses are also negligible. 

The magnetic field needs to decrease greatly outside the guide, because the laser 
cooling will not work well if the field exceeds 1 mT. It must also decrease rapidly because 
the cloud of molecules expands as soon as it leaves the guide. A steel annulus (5 mm 
inner radius and 2 mm thick) attached to the end of the guide provides a good fiux 
return path between magnets and forces the field to drop to an acceptable level only 
4 mm from the end of the permanent magnets. 

4. Laser Cooling 

Figure [2] shows a simplified version of the levels involved in laser cooling YbF, together 
with the wavelengths of the lasers required to excite the three X^S+(A^ = 1) states 
having = 0, 1, 2. In reality, each of these states has four hyperfine levels, labelled by 
((A^, 5)^, /)F, where s = | represents the electron spin, j is the intermediate quantum 

Table 1. Factors that determine the number of molecules available for the EDM 
fountain experiment. 



Number of X-state molecules/pulse in cell 


1 X 10^3 


Fraction in thermal beam 


5 X 10-3 


Fraction having {v = = 1) 


2.4 X 10-1 


Guide transmission 


6.5 X 10-2 


Fraction cooled by molasses 


7.6 X 10-3 


Fraction returning from fountain 


7.5 X 10-2 


Total number of cold molecules/pulse 


4.4 X 10^ 
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number for A^^ + 5 , / = | represents the fluorine nuclear spin, and F the total angular 
momentum. There is no Yb nuclear spin because we use the even isotope ^''^Yb. Table[2] 
shows the hyperflne energies of the v = Q and v = 1 manifolds [43j . The v = 2 hyperflne 
intervals have not been measured, but they will not differ greatly from those in v = Q 
and V = 1. Because these intervals are large compared with the 6 MHz natural width of 
the — X^S+ transition, single-frequency excitation of one hyperflne level simply 

pumps the molecules into the other levels and the light scattering quickly ceases. In 
order to prevent this, the lasers are modulated to produce sidebands that excite all 
four hyperflne levels. The hyperflne splitting of the upper state A^Iiii2{v = 0) is small 
enough [44j to neglect. While the light contains all the necessary frequencies, it still 
does not contain all polarisations and therefore the molecules optically pump into a dark 
superposition of Zeeman sublevels of the X-state. This problem is overcome by applying 
a magnetic fleld at an angle to the linear polarisation of the light, which continuously 
mixes the bright and dark states at a rate comparable with the optical pumping rate. 

We have modelled the laser cooling so that a realistic estimate can be made of the 
number of molecules in the fountain. The starting point is to determine the scattering 
rate of a molecule in the presence of the cooling light. For the moment we neglect decay 
to vibrational levels above v = 2; we shall return to this issue later. For simplicity, 
we assume that all the hyperflne levels are excited with a single common detuning 
5 and that the total intensity of the light, 12/, is equally distributed over all twelve 
frequency components (4 hyperflne levels times 3 vibrational states). To estimate how 
the population evolves in time over the 36 lower magnetic sub-levels, and the 4 upper 
sub-levels, we solve the rate equations for the 40-level molecule interacting with the 12 
laser frequencies, including an appropriate damping of population differences between 
pairs of states that are mixed by the magnetic fleld, B. This damping rate is only 
non-zero between states of the same j and F, that differ in Mf by ±1. Between these, 
the damping-rate is set to gfisB where the g-factors are -1/3, 5/6, and 1/2 for the 
(j, F) = (1/2, 1), (3/2, 1) and (3/2, 2) states respectively. The branching ratios from all 
upper levels to all lower levels are calculated using the appendix of [45j together with 

Table 2. Hyperfine energy levels of the state X^I]+(A/' = 1) for vibrational states 
V = and v = 1. 



State (i, F) v = {MHz) v = 1 (MHz) 



(3/2,1) 192.1 205.8 

(3/2,2) 161.2 176.5 

(1/2,0) 155.7 151.9 

(1/2,1) 
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the Franck- Condon factors /o_o = 0.928, /o-i = 0.069, and /o_2 = 0.0027, taken from 



Recall that the scattering rate for a 2-level system is given by ij^iii^l^l(^s/vy \^ 
where F is the spontaneous emission rate of the upper state and /gat = 7r/icF/(3A^) is 
the saturation intensity (4.4mW/cm^ for our transition). Despite the complexity of the 
real 40-level system, one finds that the scattering rate R given by the solution to the 
rate equations is exceedingly well described over a wide range of intensity and detuning 
by a function of the same simple form: 

l + //(5/3at) + 4(Vr)2l0.5 ' ^ ^ 

where / is the laser intensity in each of the 12 frequency components. The width of the 
scattering resonance is still F but the limiting rate at high intensity is reduced to F/10.5, 
rather than the F/2 of a two-level system. This is close to what one might expect - for 
strong saturation the population is distributed equally over all 40 levels, leaving 1/10 of 
the population in the A-state, rather than the 1/2 of the two-level system. The small 
difference between 10 and 10.5 is probably due to the dark states which slow down 
the scattering a little while they are rotated into bright states by the 0.3 mT magnetic 
field. The 5- fold increase in the eflFective saturation intensity similarly refiects the higher 
multiplicity of the ground state compared with the upper state. See Appendix A for 
further discussion of this. 

The simple form of ([T]) allows a straightforward extension of standard Doppler 
cooling theory [37j to our case, where the saturation parameter s = ///sat is replaced by 
5eff = //(5/sat) and F is replaced, where appropriate, by Feff = j^F. In a ID molasses, 
the coefficient of the frictional force —av is then 

(l + .eff + 4(VW r2 

and the minimum temperature. 



^=7^ . ./c.^^o^2 ' (2) 



'''''' 2kB 

is reached for a detuning of 



Vl + '^eff , (3) 



(^Tmin = - 2 Vl + '^eflf • (4) 

These formulae provide analytical insight into the operation of the molasses, but 
cannot tell us the number of ultracold molecules it produces. To determine this we 
have done a three dimensional numerical simulation of molecules moving through the 6- 
beam molasses. We use the effective excitation rate ^ i+{25/vY ^ effective decay 

rate Feff, so as to reproduce the scattering rate given by equation ([T]). The detuning 
5 includes the Doppler shift Vi/\ where Vi is the velocity component in the direction 
of beam i. The simulation tracks the position and velocity of each molecule by taking 
account of the momentum it exchanges with the light whenever a photon is absorbed 
or emitted. 
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Figure 5. Phase space plots showing the distribution of molecules in a plane transverse 
to the guide axis as they are cooled by the optical molasses. 



The molecules emerge from the lOmm-diameter guide described in Sec.|3] at a 
temperature of 3.6 K, with the transverse velocity truncated at ±5.9 m/s. The molasses 
is centred 1.3 cm downstream from the end of the permanent magnets. Each light 
beam has an intensity profile e~^^/^^ with w = 1 cm and each laser has a total power of 
275 mW. The laser power is recycled to form the 6 beams. The detuning from resonance 
is 5 = — 27r X 5 MHz. We first simulate a set of molecules that have the full distribution 
of velocities but all exit the guide at the same moment in time. Figure [5] shows the 
distribution of molecules projected onto the x — plane of phase space at several 
moments of time since they emerge from the guide (x is the vertical direction in figure 
[3]). Initially, the acceptance of the guide is filled and the phase space density is almost 
uniform. After 5 ms, the transverse velocity is strongly compressed corresponding to 
ultra-low temperature. Similar results are found in the other two directions. From the 
simulation we find that the molasses has a capture velocity of approximately 10 m/s. 

Once the molecules enter the molasses a fraction of them begin to cool to low 
temperature. Figure [6] shows how the temperature of this cold fraction decreases over 
time, as determined by fitting their velocity distribution to the Boltzmann distribution. 
We see that the temperature reduces to 185 /xK, close to the Doppler limit of 194 /xK 
given by equation ^ for the parameters used in the simulation. It takes about 7 ms 
for the molecules to reach this base temperature. It is likely that the temperature will 
actually go lower than that because of additional cooling by the Sisyphus mechanism 
[37J, which is not included in our simulation. Conservatively, we assume for the moment 
that the temperature is 185 /xK. 

We turn now to the probability of decay into the = 3 ground state. According to 
the simulation, the mean scattering rate in the molasses is about 1.5 x 10^ photons per 
second. Taking the probability of decay into = 3 as 3 x 10~^, slightly below the upper 
limit determined in [36], we see that molecules are pumped into = 3 in only 2 ms - 
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Figure 6. Temperature of the molecules versus the time since they exit the guide. 



less than the 7 ms it takes to cool them. For that reason, an additional laser returns the 
= 3 population to the cooling cycle by spontaneous decay via the v = 1 74-state, as 
indicated in figure [2j This repumping can be much faster than 2 ms, allowing the full 
scattering rate for cooling to be maintained for much longer. Eventually, however, the 
molecules will decay into = 4, so a balance needs to be struck; the = 3 repumping 
rate needs to be high enough that the molecules are adequately cooled, but not so high 
that the population is pumped into v = 4. 

The molecules we aim to capture are those with forward speeds between 2m/s and 
lOm/s, which emerge from the guide over an interval of approximately 100 ms. The 
molasses must therefore be applied for that long. We estimate that the branching ratio 
to > 3 is approximately (0.07)^ = 2.4 x 10~^, based on the trend of the Franck-Condon 
factors measured in [36]. If the molasses were fully on for 100 ms, 70% of the molecules 
captured would subsequently be lost to the higher vibrational states. This is avoided 
by repumping more slowly so that the molecules spend some time in the ^ = 3 state, 
where they are not excited by the cooling lasers. 

We have simulated this more complicated cooling situation. The molecules entering 
the guide are assumed to have a time distribution t^e"^^/^^^*^ with at = 10 ms. Their 
distribution of arrival times at the exit of the guide is determined from this distribution 
and the distribution of forward speeds. They enter the molasses where the parameters 
are the same as before, except that now a very weak ^ = 3 repump laser is added. 
For each molecule the simulation proceeds as follows. The molecule scatters n photons 
before decaying to v = 3, where n is drawn at random according to the geometric 
probability distribution rs{l — r^)^ with ra = 3 x 10""^. It remains in v = 3 until, after 
a time T, it is repumped into the cooling cycle. This time is chosen at random from 
the exponential distribution e"^/'^^, where we set the characteristic repumping time to 
Ts = 8 ms. This cycle continues until the molasses is turned off. A molecule is lost to 
'L' = 4 if it has scattered more than m photons in total, where m is drawn at random 
from the geometric probability distribution r4(l — r^)^ having r^ = 2.4 x 10~^. Of 
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Figure 7. Fraction of ultra-cold molecules versus time since the molecules were 
produced in the source. Here ultra-cold means that v < 3, that the speed is smaller 
than 30cm/s, and that the molecule is within 2 cm of the centre of the molasses. 



all the molecules that are cooled, 35% are lost in this way. There is also loss through 
spontaneous vibrational decay within the X state. This changes the rotational state 
from A^ = ltoA^ = 0or2, causing the molecule to decouple from the cooling cycle. 
The lifetimes of the excited vibrational states are given in figure [2] We calculated these 
using the dipole gradient dji/dR = 59Debye/nm [47j. We find that a further 19% of all 
the cooled molecules are lost this way. 

Figure [7| shows the fraction of molecules cooled to ultra-cold temperature as a 
function of the time since they were produced. This is given as a percentage of all 
the molecules exiting the guide. The "ultra-cold" molecules are taken as those having 
< 3, a total speed less than 30cm/s and a displacement of less than 2 cm from the 
centre of the molasses. Losses due to decay between the X state vibrational levels are 
included. Molecules begin to load into the molasses after about 30 ms when the fastest 
ones that can be captured start to arrive. This ultra-cold fraction increases to 0.76% at 
t = 100 ms, when the molasses is turned off. Holding the molasses on for longer does 
not substantially increase this fraction. The temperature of the ultra-cold distribution 
is 185 //K and the spatial distribution has an rms width of approximately 0.8 cm in each 
direction. 

Because the v = 3 repumping is so slow, the required laser intensity is low: we 
estimate ^ 30/xW/cm^ on each of the four hyperfine components. Instead of adding 
rf sidebands to address each component, a broadband laser could be used, or a single 
frequency laser could be scanned back and forth over the hyperfine structure. These 
options make the v = 3 repumping particularly simple to achieve. 
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Figure 8. Hyperfine levels of the YbF ground state X^i;+('u = 0,A^ = 0) in the 
presence of an electric field. Molecules initially in the F = state are excited by an rf 
pulse to a coherent superposition of the F = l,mi? = ±l states. 



5. Molecular Fountain 

Once the molecules are cooled, moving molasses provides a convenient way to launch 
them. The upward- and downward-going laser beams are detuned by a frequency 
±v/\ so that the rest frame of the molasses moves upward at a velocity of v (where 
the Doppler shift cancels the frequency difference). Next, the beams that excite the 
state (v^j^F) = (0, are switched off and replaced by a single retro-reflected 
beam, circularly polarised along the magnetic field. This excites the Zeeman sublevels 
rriF = —1,0 but not tuf = +1, which becomes dark. All the molecules are optically 
pumped into this dark state. All the light is then turned off, leaving the cloud to fly 
freely upward. With a velocity of 1.5 m/s, the molecules reach a height of 11 cm before 
falling down again under gravity, returning to the optical interaction region after 300 ms. 
The cloud expands as it flies, so not all the molecules return through an aperture at the 
bottom of the trajectory. With a 1cm x 4cm aperture, formed by the electrodes of the 
EDM experiment, the Maxwell-Boltzmann distribution gives a transmission of 7.5% for 
a temperature of 185 /xK and a flight time of 300 ms. This fraction of the cloud returns 
to the optical molasses region, where it is detected by laser-induced fluorescence using 
low-intensity optical molasses. Because each molecule is able to scatter thousands of 
photons, they are all detected. This constitutes the molecular fountain that will be used 
to measure the electron EDM. 

Table [T] collects together the various factors discussed above that determine the 
number of cold molecules available for the electron EDM measurement. This shows 
that the fountain can provide 4.4 x 10^ detected molecules per beam pulse in the state 



X2S+(^ = 0, = 1) with (j, F, tuf) = (|, 1, +1). 



6. Electron EDM measurement 



As the molecules fly up, a microwave vr-pulse at 14 GHz drives them to the absolute 
ground state (A = 0, F = 0), indicated in figure |8| The magnetic field is then turned off 
while the molecular cloud flies up into the electric field plates of the EDM experiment. 
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illustrated schematically in figure [3| Once it is there, we raise the field from zero to 
E ^ lO^V/m. This allows the cloud to enter the plates without being accelerated 
(to 5.5 m/s) and defocussed by the gradient of the fringe field, which would damage 
the sensitivity through the loss of molecules. The plates also form part of a TEM 
transmission line operating at 170 MHz, as demonstrated in [46j. A vr-pulse of rf power 
transfers all the molecules to the F = 1 hyperfine state. Because the rf magnetic 
field is perpendicular to which defines our z axis, the molecules are driven into 
the superposition state = +1) + \mF = —1)), as illustrated in figure 8, This 

corresponds to an alignment of the spin along the x-axis. A weak static magnetic field 
B is also applied along z. Under the magnetic and electric dipole interactions, the 
spin precesses in the xy plane at an angular frequency of {ihbB + d^rjE)/?!^ where the 
first term is the normal Zeeman interaction and the second is the small electron EDM 
interaction of interest, magnified by the enhancement factor 77. 

The molecules now fiy freely for a time r ^ 250 ms, while they rise up to 
the top of the plates and fall down again, accumulating a total precession angle 
= {i^bB + der]E)T/h. (The blackbody radiation has no appreciable effect on the 
molecules over this amount of time [48j). A second rf vr-pulse repopulates the F = 
state with probability cos^ 0. The static electric field is switched off and those molecules 
that have been returned to F = are then driven back to = 1 by a second microwave 
pulse, and detected using the optical molasses to induce fiuorescence. The fiuorescence 
signal is proportional to the number of molecules and hence to cos^ 0. A reversal of 
the electric field changes (/) to (/x^S — der]E)T/h^ leading to a small change in the 
fiuorescence intensity. This change is maximized by setting the magnetic field to 35 pT, 
so that (/) ^ tt/A (or — 7r/4 if B is reversed). 

Thus the electron EDM is measured by recording the molecular fiuorescence for 
each beam shot, together with the signs of E and B. After each shot, the apparatus 
is switched to a new configuration and the next fiuorescence pulse is measured. The 
EDM phase defjEr/h is derived from many such measurements by adding or subtracting 
the fiuorescence signals according to the sign of the product EB. In practice, it is very 
beneficial to modulate more parameters, as detailed in [49j, but this is the basic principle. 
The statistical uncertainty in d^ is 
_ 100 h 

~ ^ 2r]ETVN ' ^ ^ 

where the first factor converts from the SI units of Cm to the conventional e.cm. In 
the second factor, rjE = 1.4 x 10^^ V/m is the effective electric field interacting with the 
electron EDM when the applied field E is 10^ V/m [13j, r is the free precession time and 
is the total number of molecules detected in the experiment. With 4.4 x 10^ molecules 
per shot (Table [T]) and two shots per second, the statistical noise will be 6 x 10~^^ e.cm 
in 8 hours, the level indicated by the solid arrow in figure [T] 

In order to achieve sensitivity at this statistical limit, the noise due to random 
fiuctuations of the magnetic field must be suppressed to below 25fT/v/Hz, which can 
be achieved with good magnetic shielding, together with the use of appropriate materials 
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inside the apparatus. Excess noise due to non-statistical fluctuations in molecule number 
can be removed by normalising to the fluorescence intensity detected in the initial 
molasses. 

Systematic errors will also need to be under control. References [T5l Il9] describe all 
the signiflcant systematic errors that are known, and strategies for dealing with them. A 
leading concern will be the magnetic fleld change associated with electric fleld reversal. 
For this reason, leakage currents will be reduced to the pA level, as demonstrated in [50j. 
Even when there is no leakage, the currents that flow when the plates are charged (about 
10 /xA if the fleld is switched in 10 ms) can magnetise the shields. A design to reduce 
this adequately is discussed in [5Tj. It will be important to show by direct measurement 
that the magnetic fleld is properly behaved. For this purpose, we plan to use flbre- 
coupled SERF magnetometers [52]. This recently-developed device is ideally suited to 
our needs: it has a fleld sensitivity of a few fT/v^Hz and operates at flelds below 10 nT. 
Another systematic error considered in [15j arises from the geometric phase [53] due to 
the (slight) rotation of the electric fleld viewed from the rest frame of the molecules. In 
[15] this gave an uncertainty of 3 x 10~^° e.cm. In the proposed apparatus, the geometric 
phase will be of similar size but the coherence time will be 390 times longer, making the 
EDM uncertainty correspondingly smaller, and thus entirely negligible. 

7. Summary 

In conclusion, we have shown that laser cooling of molecules will make it possible to 
build an intense, cold molecular fountain - an important new tool for quantum metrology 
and fundamental physics studies. We have outlined how such a source may be used to 
improve experimental sensitivity to the EDM of the electron by up to three orders of 
magnitude. This speciflc application is important because it probes new elementary 
particle physics at energy scales up to 100 TeV and promises to illuminate the CP- 
violating physics of the very early universe. 
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Appendix A. Further discussion of excitation by the molasses lasers 

Consider a molecule with excited state e that is coupled by the laser beams of the 
molasses to ground states Qj. Let the probabilities of occupation, rig for the excited 
state and rij for the ground states, be governed by the rate equations 

d n ' 

= AjTie + Rj{ne-nj), ( A. 1 ) 
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where Aj is the partial spontaneous decay rate from e to gj and Rj is the excitation 
rate for that transition, together with the normahsation condition 

N 

He + = 1. (A. 2) 

In the steady state, where dnj/dt = 0, these equations give an excitation probabihty of 

rie = \ . (A.3) 

{N + 1) + EA,/R, 

The ratio Aj/ Rj can be written as 2Isat,j/Ij^ where Ij is the intensity of the hght driving 
the j^^ transition. Igatj = TrhcT / {3X^) is the saturation intensity for that transition, Aj 
being the wavelength. F is the total decay rate of state e. This gives 

ne = . ( A.4) 

As the wavelengths involved in the laser cooling do not vary greatly from one vibrational 
state to another, it is instructive to approximate the Igatj as a constant Igatj which can 
be brought outside the sum to give 

Tie ^ . (A.5) 

(N + l) + 2Isat E 1 A. 

In order to have the strongest cooling, one wants to make rig, and hence the scattering 
rate, as high as possible given the available laser intensity. It is most efficient to choose 
equal values for all the /j, making Ij = Itot/N ^ where Itot is the total intensity of the 
light. Then the excited state probability becomes simply 

^ (iV + 1) + 2N^Isat/Itot ' ^^-^^ 

and the ground state populations rtj are all equal to (1 — ne)/A^. When Itot « j^hat^ 
the intensity is low and tends to {l/N'^)Itot/{'^Isat)- At the opposite extreme, when 
the intensity is high, the population is shared equally among all the levels and tends 
to 1/(A^ + 1). 

In the case of YbF there are 36 lower states, these being the twelve sublevels of 
F = 2, 1, 1, multiplied by the 3 vibrational states v = 0, 1, 2. Instead of one upper 
level, there are four: the sublevels of F = 0, 1. At any given moment, each excited state 
is only coupled to a subset of the ground states, those having a non- vanishing transition 
matrix element. However, we apply a magnetic field that mixes light and dark states so 
that every state can be excited. Our numerical simulation yields an empirical form for 
the scattering rate, given in equation ([T]). After dividing by spontaneous emission rate 
F, setting detuning 5 to zero, and writing Itot = 12/, this equation gives the empirical 
excited state probability as 

4 

42 + 2 X 35.52/,,,//,,, • 
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Comparison with equation (A. 6) shows that, despite the complexity of the real system, 
it is not very different in the end from four upper levels, each separately connected to 
36 lower levels. 
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